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We present a method for the reahzation of radial and azimuthal polarization for non-paraxial beams in a 
rigorous but simple manner. This is done by constructing independent exact vector solutions of Maxwell's 
equations from scalar Bessel beams and combining them together by analogy with the paraxial case. In this 
way we are able to obtain spatial and polarization patterns analogous to the traditional cylindrically polarized 
paraxial wave field. Applications of these beams are discussed. © 2013 Optical Society of America 
OCIS codes: 260.0260, 260.5430, 260.6042 



C/5 

o 

• I— I 

> . 

Oh- 

q 
o 

• 1— I 
Oh- 



> 

On 
(N 

O 

m 



X 



Among the broad zoology of solutions of Maxwell's 
equations, radially and azimuthally polarized vector 
beams, i.e., beams of light with non uniform polarization 
patterns, attracted great interest in the last decade [1-3]. 
Their particular character, in fact, makes them very ap- 
pealing for applications in various fields of research such 
as single molecule spectroscopy [4], optical tweezing [5], 
confocal microscopy [6] as well as in commercial applica- 
tions like material processing [7,8] and data storage [9], 
just to name a few. Another interesting feature of such 
vector beams is that they can yield to very small focal 
spots [10] and generate axial electric [11] or magnetic [12] 
fields once focused. All these reasons have motivated the 
development of several experimental techniques to real- 
ize such field configurations [13-16]. Recently, it has been 
also shown that squeezed azimuthally polarized optical 
beams show hybrid entanglement [17]. For a comprehen- 
sive theoretical analysis of the properties (both classical 
and quantum) of general cylindrically polarized states of 
light, the reader is addressed to Ref. [18]. 

Despite the considerable amount of work that has been 
done in this field, radially and azimuthally polarized 
beams have been studied only within the framework of 
paraxial optics. The only exception, to the knowledge of 
the authors, is represented by the work of Bouchal and 
Olivik [19], where the concept of non-diffractive vector 
beam is introduced and extended to the non-paraxial do- 
main. However, we feel that the physical meaning of the 
solutions found in Ref. [19] tends to be obscured by the 
use of a heavy mathematical formalism. It is our aim 
here to introduce cylindrically polarized states of light 
beyond the paraxial domain in a much simple way, by 
analogy with the paraxial case which leads to a clearer 
physical understanding, thus making this important sub- 
ject accessible to a broader audience. Thanks to the fact 
that we use Bessel beams as basis to create these new 
fields, they will also possess the typical characteristics 
of Bessel beams, such as self healing and diffractionless 



propagation. They could also be used, as it has been re- 
cently proposed [20], to generate accelerating beams of 
arbitrary trajectories. All these features, combined with 
the complex cylindrical polarization pattern introduced 
here, may lead to new possible developments in timely 
research fields such as strongly focused light or to inves- 
tigate the role of complex polarization patterns in light- 
atom scattering problems [21]. 

Let us start our analysis by considering a monochro- 
matic scalar field '^(r) at angular frequency uj that is a 
solution of the scalar wave equation (V^ + k'^)ilj{T) = 0, 
where k = cj/c is the wavenumber. It is not difficult 
to demonstrate [22] that the time-harmonic electric and 
magnetic vector fields are expressible in terms of the 
scalar field ip{r) as 



E(r) iuj 



V(V • A) 



B(r) 



V X A, 



(la) 
(lb) 



where A = £^o'0(r)f , being f = ax + /3y, with |ap + 
= 1, is a complex unit vector orthogonal to the 
propagation axis z that defines the polarization of the 
beam and Eq is arbitrary. For the sake of clarity, in the 
remaining part of the manuscript, we will focus our at- 
tention on the electric field solely. 

As shown in Ref. [18], radially and azimuthally polar- 
ized beams live in a four dimensional space spanned by 
the basis formed by the Cartesian product of the base of 
Hermite-Gauss modes ?/^nm(r) of order = n + m = 1 
[23] {V^io, V^oi} and the polarization basis {x, y}, namely 
{V^io, V^oi} {x,y} = {V^iox, V^ioy, V^oix,?Aoiy}. Linear 
combinations of these four vectors give rise to radially 
{ur) and azimuthally (ua) polarized fields as follows: 
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^(^V^iox + V^oiy), 



^(tV^oix + V^ioy), 



(2a) 
(2b) 



1 





Fig. 1. Contour plot (a) and three dimensional intensity 
profile (b) of the scalar function ^io(r). Panel (a) shows 
how, close to the propagation axis, ^io(r) correctly re- 
produces the behavior of the Her mite- Gauss beam i/jio 
employed in the paraxial case. Note from panel (b) that 
for x/ Wo > 4 and y/wo > 4, the scalar function ^io(r) 
shows some ripples, whose intensity is much smaller than 
the central lobes. 



where the ± sign indicates co-rotating and counter- 
rotating modes respectively [18]. It is well known in 
paraxial optics that suitable linear combinations of 
Hermite-Gaussian modes may lead to orbital angular 
momentum carrying beams [24]. In particular, it is pos- 
sible to demonstrate that the Hermite-Gauss modes 
V^io and t/jQi can be written as a superposition of two 
Laguerre- Gaussian beams with angular momentum m = 
±1 [23]. Having this in mind, we now construct the 
non-par axial counterpart of the Hermite-Gaussian ba- 
sis modes, using combinations of vector Bessel beams in 
such a way that the angular momentum of the resulting 
non-paraxial beam is preserved. 

First, we choose for the scalar field a 
mth-order Bessel beam, namely i^mi^) — 
Jm{KoR) exp {imcj)) exp {iz^^/k'^ — Kq), where 
Kq = ksinOo (being the cone angle aperture of 
the Bessel beam), x = Rcoscj) and y = i^sin^. Then we 
construct, by analogy with the paraxial case, the non 
paraxial scalar beams ^io(r) and ^oi(r) by combining 
two Bessel beams ^^^(r) with angular momentum 
m = ±1 respectively, namely ^io(r) = (V^f + V^:?]^)/a/2 
and ^oi(r) = -^(V^f - V^:?i)/\/2. In Fig. 1 the intensity 
profile of the scalar beam of ^io(r) is plotted as a 
function of the dimensionless coordinates x/wq and 
y/wo being wq = 2/{k6o) the beam waist. From this 
figure it is clear that, if we limit ourselves to the central 
domain (Fig. 1 (a)), the intensity distribution of the non 
paraxial beam ^io(i*) is fully equivalent to the intensity 
distribution of the correspondent Hermite-Gauss beam 
i/jio. Discrepancies start to arise in the non-centered 
region (Fig. 1 (b)), where the Hermite-Gauss beams 
are zero because of the Gaussian damping, while the 
beam ^io(r) presents ripples that have the symmetry 
structure of the central lobes, but a sensibly smaller 



Fig. 2. Complex polarization patterns of (a) and 
(b) U^, with superimposed the donut-shaped intensity 
distribution. The axes of both graphs span the interval 
[—5, 5] in units of the beams waist wq. 



intensity. The presence of these ripples is a signature of 
the non-paraxial character of these beams. As last step, 
following Eqs. (2) we build the non-paraxial radially 
and azimuthally basis as 

U± = -i=(±*iox+*oiy), (3a) 

U± = -^(T*oix+*ioy). (3b) 

These modes are exact vector solutions of the full 
Helmholtz equation, since the scalar fields ^io(r) and 
^01 (r) are solutions of the full Helmholtz equation, in 
contrast with the case of Gaussian beams that are so- 
lution of the paraxial equation solely. However, they 
are not exact solutions of Maxwell's equations, as they 
do not evidently satisfy the transversality condition 
V • A ~ ^' Exact solutions of the Maxwell's equa- 
tions are obtained by substituting and into 
the right side of Eq. (1), to generate electric and mag- 
netic fields that are automatically satisfying the condi- 
tion V • E(r) = 0. 

The polarization pattern as well as the intensity dis- 
tributions of these non-paraxial radially and azimuthally 
polarized beams are displayed in Figs. 2 and 3. Although 
they retain their typical donut-like intensity distribution, 
a closer inspection of the polarization pattern of these 
beams reveals that they possess a nonzero longitudinal 
field component, as can be understood by noting that, 
in contrast to their paraxial counterpart (see Ref. [18]) 
they manifest a nonzero polarization in correspondence 
with the center of the intensity distribution. The inten- 
sity distributions of these non paraxial vector beams are 
displayed in Fig. 4 for both the cases of radially (upper 
panels) and azimuthally (lower panels) polarized beams. 
In particular, the panels (b) and (d) show the inten- 
sity distribution of the longitudinal component of such 
beams. As can be seen, for the radially polarized beam 
the z-component of the electric field different from zero, 
although about three orders of magnitude less intense 
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Fig. 3. Complex polarization patterns of (a) and 
(b) U^,witli superimposed the donut-shaped intensity 
distribution. The axes of both graphs span the interval 
[—5, 5] in units of the beams waist wq. 



than the transversal intensity distribution. On the other 
hand, the z-component of the electric field for the az- 
imuthally polarized beam (panel (d)) remains exactly 
zero. 

In conclusion, we have shown how to generalize the 
theoretical framework of radially and azimuthally po- 
larized beams to the case of non-paraxial fields, using 
Bessel beams as generating beams. We shown that these 
non-paraxial fields retain the typical donut-shape of their 
paraxial counterparts, although they show a nonzero lon- 
gitudinal electric field component for the case of radially 
polarized beams. 



Fig. 4. Intensity distributions for the radially (upper 
row) and azimuthally (lower row) polarized non parax- 
ial electromagnetic fields. Panels (a) and (c) display 
the transversal components of the field (|£^a:p + l^^yP)? 
while in panels (b) and (d) the longitudinal component 
X 10^ is depicted. The coordinates are expressed 
in units of wq. The intensity of the 2:-component of the 
electric field has been rescaled to the maximum of the 
transverse intensity distribution. 
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